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Abstract. Positronium (the bound state of electron and positron) has been thought 
to be formed after proton decay (> 10 34 yr) through collisional recombination and 
then decays by pair annihilation, thereby changing the matter content of the universe. 
We revisit the issue of the formation of positronium in the long-term future of the 
universe in light of recent indication that the universe is dominated by dark energy 
and dark matter. We find that if the equation of state of dark energy w is less 
than —1/3 (including the cosmological constant w — —1), then the formation of 
positronium would not be possible, while it is possible through bound-bound transitions 
for — 1/3 < w < —0.2, or through collisional recombination for w > —0.2. The 
radiation from e^- pair annihilation cannot dominate over e , while that from proton 
decay will dominate over baryon and e for a while but not over dark matter. 
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1. Introduction 

Since the pioneering work by Dyson |T] (see also |2j), the physical processes that would 
occur in the far future of the universe have been studied by several authors [3 BJ El Hj ■ 
One of interest then was the effect of proton decay on the matter content of the universe. 
It was found for an open universe that the decay of protons (with lifetime td) will produce 
nonthermal radiation which will dominate the universe from tj, to lOt^ and then matter 
domination follows Moreover, electrons and positrons from proton decay will bound 
to become positroniums and eventually annihilate [3J. However, the annihilation of 
positron-electron pair occurs so slowly that the radiation from pair annihilation cannot 
dominate the universe [3]. 

Today the situation in cosmology becomes much changed. The new ingredients 
after Page and McKee's paper are: (1) Cold dark matter; baryons are not dominant 
component of the universe, (2) Dark energy; the universe is presently dominated by 
energy density with negative pressure. As an additional minor point, we note that 
the lower-bound of the proton decay, > 10 33 yr [H], is much larger than the previous 
bound (t d > 10 30 yr)- The purpose of this paper is to revisit the nature of the long-term 
future of the universe taking into account these new ingredients after Page and McKee's 
study jl]. Throughout this paper, we assume that the laws of physics continue to hold 
unchanged. In particular, we assume that the charge of electron does not vary with 
time. 

2. Positronium Formation in an Expanding Universe 

According to the grand unified theory (GUT) of elementary forces, all baryons are 
unstable. The life time of a baryon td is predicted to be typically some 10 34 years for 
(supersymmetric) GUT models. The positrons and electrons, when first produced from 
nucleon decay, will have energies of several hundred MeV so that they will appear as 
relativistic particles until the universe has expanded by a factor of the decay energy 
divided by the electron mass. Let g_ be the electron mass divided by the average decay 
energy per the electron mass and q + be the same ratio for positrons. Let /_ be the 
average number of electrons produced per nucleon decay and /+ be the same quantity 
for positrons. For SU(5) GUT model, \ assuming that 16% of the nucleons in the universe 
are neutrons, they are given by |U| 

g_ = 5.9 x HT 3 , /_ = 0.54 (1) 

g + = 2.1xl0- 3 , /+ = 1-38. (2) 

% We adopt non-supersymmetric SU(5) GUT as an example only for simplicity. Our qualitative results 
will not be sensitive to the detail of the model. 
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Thus, baryon density ps, e ±, s density p e and radiation density p r evolve according to 
§0 

p B = - 3Hp B - t^pB (3) 

p e = - 3Hp e + ^f- (q^f-6(t_ - t ) + q + f+6{t + - t )) (4) 
m p t d 

p r = - AHp r + tfpB ~ ^ {q-f-6(t. - t ) + q + f + 9(t + - t )) , (5) 

mpt d 

H 2 = -^-{pDE + PDM + PB + Pe] Pr) (6) 

where m is the mass of electron(positron), m p is the proton mass. £-(£+) is the time 
when the universe expand l/g_(l/g + ) times from now (to), and pde (pdm) is the energy 
density of dark energy (matter). The second term in Eq.(jl|) includes the effect that 
electrons (and positrons) from proton decay are initially considered radiation and will 
become matter when the universe is (l/<?+) times expanded. Assuming that the 
universe is dominated by dark energy with the equation-of-state w and dark matter with 
Pdm(^o) — 6xp B (t ) which is stable || and taking p r = and p e = (m/m p )p B (to)(f + — /_) 
(from charge neutrality) at the present time t = to, we solve the above equations. The 
results are shown in Fig. 1 for w = —0.3 and w — — 1. We find that radiation can 
dominate over baryon, electron and positron for only a short period after proton decay 
for w = —0.3. 



2.1. Necessary Condition 

When a baryon with mass m p decays, a significant fraction of its energy will go into e ± s 
(their mass m) with average rms momentum and energy p = jm at that time tj, where 
7m is roughly the fraction of energy m p of each decay that goes into each e produced, 
say 7 ~ m p /2m ~ 10 3 in accord with the previous example (l/g_,l/g + ). As the 
universe expands, in the absence of non-gravitational interactions the rms momentum 
will redshift to jm(ad/a). Once this becomes less than m (when a > "fad), the average 
kinetic energy per e will be p 2 /2m ~ / ~f 2 m(ad/a) 2 . In the following, we shall consider 
a cosmological model with dark energy with equation-of-state w. Hence the expansion 
law is a oc t 2 ^ 1+w \ 

As Barrow and Tipler argued, a positronium will form if the average energy of an 
electron E at average distance r from a positron 

E ~ "f 2 m(a d /a) 2 - e 2 /r (7) 

§ Here we have assumed that all of the baryons in the universe are free at the epoch of proton decay. 
In reality, however, a large fraction of the baryons would be locked up in degenerate remnants (mostly 
in white dwarfs) at this future time. When protons decay within dwarfs, the positrons annihilate with 
electrons, and the resulting high energy photons are reprocessed into low energy photons. As a result, 
protons that are locked up in degenerate remnants are essentially lost to the processes discussed here. 
The subsequent analyses would apply to the baryons that will not be processed into stars and stellar 
remnants. ^3 13 

I We assume that dark matter is stable. The subsequent discussion will be greatly changed for unstable 
dark matter, which will be strongly model-dependent however. 
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Figure 1. Evolution of ps, p e and p r for w = —0.3 (solid lines) and for w = — l(dashed 
lines). a corresponds to the scale factor at the present time t a . p e does slightly increase 
after baryon decay although it is barely seen in the figure. 



is negative 3\.% Note that the Coulomb energy decays slower than the kinetic energy 
and hence E will be eventually negative. Here the mean separation r is estimated from 
the number density of e ± , N 

r ~ iV~ 1/3 . (8) 

iV is related to the number density of baryon at the decay time Nd as 

N ~ N d (a d /af. (9) 

% This is the case for bound-bound transition. Positronium can be formed much earlier through 
collisional recombination 0|. The detailed anaylsis is given in Appendix. 
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Nd is written in terms of the ratio of baryon number density to entropy density 
n B /s ~ 10- l0 {Q B h 2 /10- 2 ) as 

N d ~ n B0 (a /a d ) 3 ~ 10~ 7 (r^s" 1 /^ 10 ) (a /a d ) 3 cm- 3 , (10) 

where a is the scale factor at present. An electron-positron pair will bound if E < 
when a > a*, where from Eq. (J7J) 

a*/a d ~ frne-V' ~ lO^lO^^^/lOV) 273 ^- (H) 
In an accerelating universe (w < —1/3), another condition for the formation of 
positronium arises: repulsive gravitational force sourced by dark energy should not 
overcome the attractive Coulomb force: 

e 2 /r 2 > —mra/a = — mr(l + 3w)H 2 /6. (12) 

The question is whether Eq. (|12J) is satisfied at that time. From 

r ~e 2 7 - 2 ^ _1 K/^) 2 (13) 

and 

H- 1 ~ H \a/a f^ 2 ~ 10 34 (t d /10 34 yr)(a/^) 3(1+, " )/2 yr, (14) 
it is required that (for w < —1/3) 

aja d < lO^-^/lOV) 273 ^" , (15) 

where we have neglected an unimportant numerical coefficient coming from |1 + 3w\ in 
Eq.(|12jl. Two conditions Eq. lflljl and Eq. p5jl are compatible if (for w < —1/3) 

10 21 10 16/(l+ W ) (t(j/10 34 yr) 2/3(l+«) < 10 49/(l-«,)( f<j/10 34 yr) 2/3(l-«,)_ (1 q) 

Taking the logarithm, the above condition is rewritten as 
4w , , 34 s 33 + 65m; 



log(t d /10 3 V) + -; r - 21 > 0. (17) 

3(1 — w z ) 1 — w z 

It is easily found that the left-hand-side of Eq. (jl7j) is an increasing function of w (for 
— 1 < '«;<— 1/3) and negative as long as 

w < -1/3 (18) 

for td > 10 33 yr + Hence positronium will not form in the universe dominated by dark 
energy with w < —1/3.* We will restrict ourselves to the case with w > —1/3 hereafter. 

The formation of positronium through collisional recombination is studied in 
appendix. Here we only give the results: the typical time scale for the formation is 
~ 10 83 yr(^/10 34 yr) and the size of positronium is ~ 10 45 Mpc for w = —0.1 (see Fig.3 
and Fig. 4). 

+ The condition would apply to the equation of state at the time of formation. For dark energy which 
has the present equation of state w < —1/3 temporarily and then would increase toward w > —1/3 
(for example, quintessence axion the formation of positronium would be possible. 

* This does not mean that positronium will never form in the universe with w < —1/3. Although some 
small bit of activity always takes place, it will be exponentially suppressed. |12j 



Ps and Dark Energy 



6 



3. Positronium Decay and Annihilation 

We have seen that positronium will be formed if w > —1/3. Positronium, once formed 
long after nucleon decay, will decay toward ground state and then pair-annihilate. One 
can use the corresponding principle to calculate the time within which the positronium 
will spiral into the ground state by dipole radiation (for large n) [THl Ej 

tspirai ~ m~V 10 n 6 . (19) 

From Eq. (|f)2jl in Appendix up to n max , the fraction of free e^s evolves as 

y e oc a~ 3(2 - Q)/5Q , (20) 

where the expansion law is a oc t a , a = 2/3(1 + w). Then a lifetime of decay to the 
ground state is 

t sp irai oc n 6 max oc y^a 3 oc a 6 ^ 1+2a ^ 5a . (21) 
Expressing y e in terms of t sp i ra i gives the logarithmic annihilation rate coefficient |4. 

t 2 + Aat 7 + 3wt' ^ ' 

Thus, after the formation of positronium, e ±, s density p e and radiation density p r evolve 
according to 

p e = - 3H Pe - pr l p e (23) 
p r = - AHp r + f3r 1 p e . (24) 

(25) 

Assuming that p e = p ea - and p r = p r d at t = td, the solutions are 

Pe = Ped{t/t d )-^ + V (26) 
Pr = ^ {(t/t d )-^ - (t/t d )" 4a ) + p Td {t/t d )- Aa . (27) 

a — p 

The asymptotic ratio of p e to p r is given by 

Pe _ a- (3 __ 8 - 9w - 9w 2 

p r ~ (3 ~ 3(2 + 3w)(l + w)' ^ ' 

In Fig. 2, the ratio is shown as a function of w. The ratio is greater than 4/3 for w < 0. 
Hence, the radiation by pair annihilation cannot dominate over e ± . This is because 
in the universe with dark energy the expansion rate is higher and radiation suffers from 
much larger redshift. 



4. Summary 

We have investigated the matter content in the long-term future of the universe. After 
the proton decay, radiation becomes dominant over baryon and e ± for a while and then 
is followed by them again but it will not dominate over dark matter. 
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now Dark Energy Dark Matter Baryon Radiation e 

10 34 yr Dark Energy Dark Matter Radiation Baryon e 
10 38 yr Dark Energy Dark Matter e ± Radiation 
10 83 yr Dark Energy Dark Matter e ± /Ps Radiation 



Table 1. The energy content of the universe assuming t d = 10 3 years. PS denotes 
positronium. 
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We have also investigated the possibility of formation of positronium after proton 
decay. We have found that positronium will not be formed in the universe if the equation 
of state of dark energy w is less than —1/3. Positronium will be formed after 10 83 years 
if the equation of state of dark energy is w > —1/3, and radiation will be produced 
by positronium decay. However, it will never dominate over matter due to the slow 
annihilation rate. The results are summerized in Table 1. 
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Appendix: Collisional Recombination 

In this appendix, following Page and McKee [I], we give the detailed analysis of 
positronium formation by collisional recombination. 

As shown by Page and McKee j3j, radiative processes are unimportant for free e . 
Therefore positronium will form by collisional recombination, 

e" + e + + e ± ^Ps n + e ± , (29) 

where Ps n denotes positronium with principal quantum number n. Once formed, 
positronium will tend to decay toward the ground state or at least an S state and 
then annihilate. At first, the transitions will be induced mainly by collisions (for large 
n states), 

Ps n + e ± -> Ps n / + e ± (30) 
Ps n + Ps n » -> Ps n / + Ps n /«. (31) 

Let N e and N n be the number densities of free e ±, s and Ps n 's. We introduce the 
fractions of free e ±, s and Ps n 's 

y e = N e /N } y n = N n /N, (32) 

where N is the total number density of all e s. 

The transitions Eq.([29j) and Eq.()30|) and their inverses occur at rates proportional 
to the densities of the participating particles, The rate is given by 

-J^ = ^N 2 ylC i)n - Ny e C nii y n + Ny e ^ ( a ( n '> n )Vw ~ n')y n ) , (33) 

n' 

where Ci >n is the collisional recombination rate coefficient to the nth level, C n ^ is the 
collisional ionization rate coefficient from the nth level, and a(n,n') is the collisional 
transition rate coefficient. They are defined by 

C n , = (a ion (n)v) (34) 
a(n,n) = (cr(n,n)v), (35) 
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where a ton (n) is the cross section for ionization from the nth level, a(n,n') is the cross 
section for transitions from the nth to the n'th level, and the averages are taken over 
the Maxwellian distribution of the relative velocities. By detailed balance is given 
by 

C t , n = 8n 2 (n/mTf 2 exp(me 4 /4n 2 T)C n>i , (36) 

where T is the temperature, T ~ r y 2 m{a < i/ a) 2 . Note that the binding energy of 
positronium is me 4 /4n 2 T. 

The ionization cross section of an atom is given in jTS| 

ion 7re 4 (5E + 2I n )(E - I n ) 

3EH n {E + 3I n ) ' {6I) 
where E is the energy of the incident electron and I n is the ionization energy of the atom 
(J n = me 4 /4n 2 T for positronium). However, for our crude purposes, a rough expression 
at high and low energies is adequate: 

ion £ 4 (E — I n ) . . 

a ~ WT n ■ (38) 

In terms of the ratio of ionization energy to temperature x n = I n /T, C n> i and Cj jTl are 
given by using Eq.fpHJ) and Eq. (|r?Hjl 

C n , ~ n 2 m-V 2 T-V 2 (l + x^" 1 exp(-x n ), (39) 
C itn ~n i m- 3 T- 2 (l + x n )- 1 . (40) 

The collisional cross sections for transitions between highly excited levels of 
hydrogen have been considered in several papers [T^lTi] (see [3] for a list of references). 
An interpolation formula in [Hj) gives the cross section which is accurate to 15% for 
I n < E <C m,e 2 n~ l < v <C 1. If we drop numerical factors and logarithms, then for 
1 < s e n' - n < n the cross section become 

cr(n, n + s) ~ n 4 s~ 3 m _1 i? _1 , me 4 n~ 1 s~ 1 < E < m (41) 
<j(n,n + s) ~ n 6 s _1 e _8 m~ 3 i?, me 4 n~ 2 < E < me n~ x s . (42) 

From Eq. (|35jl . we obtain 

a(n, n + s) ~ n 4 s~ 3 m" 3/2 T~ 1/2 , meV 1 ^ 1 < E < m (43) 
a(n, n + s) ~ n 6 s _1 e" 8 m" 7/2 T 3/2 , me 4 n" 2 < E < me 4 n _1 s _1 . (44) 

By detailed balance, the deexcitation rate coefficients are give by 

cn(rj + s, n) = (1 + s/n)~ 2 exp(x n — x n+s )a(n, n + s). (45) 

It is apparent that from Eq. ()43j) and Eq. ()44j) that for T > me 4 n" 2 the dominant 
bound-bound transitions are for relatively small changes s in the level n. Hence if one 
views 

a s = a s (n, t) = a(n, n + s) (46) 
y = y(n,t)=y n (47) 
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as smooth functions of n as well as of t, uses the detailed balanced relation Eq. (J45|) . and 
expands Eq. (JH3j) to second order in s = n' — n, one gets a diffusion equation 

^ ~ ^N 2 y 3 e Q,n ~ Ny e C n>i y 

d ( n d \ 

+ I Ny e ^a s s 2 n 2 exp(x n )— (exp(-x n )yn~ 2 ) J . (48) 

Neglecting a Inn factor, the sum over s gives 

n 

~ rfm-WT-V 2 , me 4 n- 2 < T < m. (49) 

8=1 

Now if we define the fraction of e s in each individual positronium state as 

z = n- 2 y(n,t)/2, (50) 

then for nonrelativistic temperatures well above the ionization temperature (T <C m 
and x n <C 1), 

dz Ny e \ 2 f Ny 2 (51) 



dt m^T 1 / 2 



m 3 / 2 T 3 / 2 J ri 1 On \ dn J 
When y e ~ 1, the population of positronium states with i„ < 1 increases toward 
the limiting value 

z s (n,t) ~ Ny 2 e (mT)-V 2 (52) 
~ N d y 2 e ( im )- 3 ~ 10- 48 1 - 48 /( 1+w )y e 2 ( td/10 34 yr) -2/(i + «,) ; (53) 

which is the Saha equilibrium value for x B < 1. The population reaches this nearly 
stationary equilibrium only for levels in which the diffusion time is short compared with 
the expansion time t. Populations are only filled up to z$ for 

n > n eq (t) ~ y^N-^^mTf'H- 1 ' 2 (54) 

~ io- 10 io 24 /( 1+w )^ 1 / 2 (t,/ioV) (1 ^ )/{2+2 " ) (V^) (1_3u ' )/(6 t^) 

Since positronium must have a orbital size r n (= 2m~ 1 e~ 2 n 2 ) less than the mean 
separation [y e N)~ 1 ^ between free e^s, 

n < n max {t) ~ m 1/2 e(i/ e iV)- 1/6 (56) 

~ 10 6 10 8 /( 1+W )y e - 1 /6( td/l0 34 yr) l/3(l + «,) ( ^ d) l/3(l +to )_ (57) 

As the universe expands, n max (t) increases and more Ps n levels become available to be 
filled up to population z$- The filling up of levels with x n < 1 decreases the fraction of 
free e ± 's by 

dy e ~ -z s n 2 max dn max ~ -yl /2 a 1,2 da/ 'a 3 / 2 \ (58) 
where the time scale a 2 is 

a 2 /a d ~ 7 2 me" 2 iV- 1/3 ~ 10 21 10 16 /( 1+ ^(Vl0 3 V) 2/3(1+w) - (59) 

This is the same as a* given by Eq.flTTJ since the binding energy, J n , is smaller than the 
average kinetic energy. 
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However, for level x n > 1 or 

n < n thermal = (me 4 /AT) 1 ' 2 (60) 
~ eWo*) ~ l^(t/t d f^ +w \ (61) 

the exponential increase in the Saha equilibrium value above zs (for small x n ) allows 
the value of z to exceed z$ if the collisional recombination rates are fast enough. Bound- 
bound transitions are less important in this case, and ionization becomes exponentially 
damped for x n > 1, so 

dz^N^yl n 4 N 2 yl n 4 N 2 yl 

m " ~8^"~ ~ {a/ad) ■ (62) 

This implies that once the ionization becomes unimportant for levels n, that level will 
fill up to population 

z ~ n 4 N 2 y 3 el ~ 2 m- 5 e-\a/a d yH. (63) 
This in turn depletes the free e ± 's by 

dye ~ -zn 2 thermal dn thermal ~ y y+^/ 2 da/a\^ 1+w)/ \ (64) 
with the time scale a±, where 

a,ja d ~ [ e - 10 7 9 m 5 iV- 2 iJ (a /a,) 3 ( 1+ -)/ 2 ] 2/{9+3w) (65) 

^ 1 Q(296+104«;)/(9+3«;)/(l+«;)^^ 1 Q34 yr ^2(3-«;)/(9+3«;)/(l+«;)_ (gg) 

In Fig. 3, t\ and t 2 are plotted as a function of w, where £i(t 2 ) is the time corresponding 
to ai (a 2 ). Typically t± ~ 10 83 yr(t d /10 34 yr)- 

For —0.2 < w < 0, ti < t 2 , Most e ± 's will bind around the time t 1: going into Ps n 
levels n somewhat below 

ni = n th ermai{ti) ~ I0' 5 (a 1 /a d ). (67) 

These positronium states are very loosely bound, with radii 

f n ~ n\e~ 2 m~ 4 (68) 
~ 10~ 43 10 (592+208w)/{9+3 "' )/(1+ "' ) (t d /10 3 V) 4(3 ~ w)/(9+3 "' )/(1+ "' ) Mpc. (69) 

r n ~ 10 30 Mpc for iy = -0.1. 

For u> < —0.2, however, since the expansion rate of the universe becomes higher, 
collisional recombination is no longer effective, and hence t\ > t 2 . The typical size of 
positronium is 

r ~ e 2 1 - 2 m-\a 2 /a d ) 2 (70) 
~ 10- 1 10 32 /( 1 - Hu )(* d /10 34 yr) 4/3(1+,o) Mpc. (71) 
r ~ 10 45 Mpc for w = —0.3. The results are shown in Fig. 4. 
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Figure 3. The formation times by collisional recombination (t\: solid lines) and 
by bound-bound transitions (t 2 : dashed lines) as functions of w. Three line are for 
td = 10 33 yr, 10 34 yr, 10 35 yr from top to bottom. 
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